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Abstract 

There are some indications from recent determinations of the strong coupling constant 
a s and the gluon condensate that the Operator Product Expansion may not be accurate 
enough to describe non-perturbative effects in hadronic tau decays. This breakdown of the 
Operator Product Expansion is usually referred to as being due to "Duality Violations." 
With the help of a physically motivated model, we investigate these duality violations. 
Based on this model, we argue how they may introduce a non-negligible systematic error 
in the current analysis, which employs finite-energy sum rules with pinched weights. In 
particular, this systematic effect might affect the precision determination of a s from tau 
decays. With a view to a possible future application to real data, we present an alternative 
method for determining the OPE coefficients that might help estimating, and possibly 
even reducing, this systematic error. 



1 Introduction 



Ever since the seminal work of Ref. pQ, it has been recognized that hadronic tau decays 
provide an excellent experimental setting for extracting precise values of Standard Model 
parameters such as the strong coupling constant a s [2], the strange quark mass m s [3], and 
the CKM matrix element \V US \ [I]. 

The importance of a precise determination of QCD parameters can hardly be overempha- 
sized, both in order to test the validity of the Standard Model as well as to constrain physics 
beyond the Standard Model. In particular, a s has a special status as the coupling constant 
of QCD. The current world average value as given by Ref. [5] is 

a 3 {M 2 z ) = 0.1189 ± 0.0010 . (1.1) 

Methods to determine a s include, among others, deep inelastic scattering, the e + e~ cross 
section, the Z width, and the lattice. What is special about the determination from tau 
decays is that it extracts a s at a remarkably low energy scale (m T = 1.777 GeV), but still 
high enough for perturbation theory to be applicable and non-perturbative effects to be small. 
Because of asymptotic freedom, the relative error on a s gets reduced by roughly a factor of 
three after evolving it from the tau mass to the Z mass. Recently, Ref. [3] performed an 
analysis based on the ALEPH data collected at LEP, quoting as the final value 

a s (m 2 T ) = 0.345 ± 0.004 exp ± 0.009 th , (1.2) 

which, extrapolated to the Z mass, gives 

a s (M 2 ) = 0.1215 ±0.0004 exp ± 0.0010 th ±0.0005 evo i (1.3) 
= 0.1215 ± 0.0012 . 

A somewhat more conservative estimate was given in Ref. [2]: 

a s {m 2 T ) = 0.340 ± 0.005 exp ± 0.014 th , (1.4) 
a s (M 2 ) = 0.1209 ± 0.0006 cxp ± 0.0016 th ± 0.0005 evo i 
= 0.1209 ± 0.0018 . 

These results are compatible with the OPAL value of a s (M 2 ) = 0.1219 ± 0.0010 cxp ± 
0.0017 th [6] . 

The results in (jl.3p . together with the lattice result [7] 

a s (M 2 ) = 0.1170 ± 0.0012 , (1.5) 

are the most precise determinations of the strong coupling constant to date. Remarkably, 
they differ from each other by 2.7 standard deviations. In fact, as emphasized by Bethke |5], 
these two values extracted from tau decay and the lattice differ from the world average com- 
puted with all the other measurements by +1.8 and —1.6 standard deviations, respectively. 
Given the relevance of these measurements, we think it is timely to assess the validity of the 
theoretical assumptions that go into the estimates of systematic errors in both cases. Our 
aim in this paper is to do this for the result extracted from tau data. 

We may subdivide the systematic uncertainties in the tau determination of the strong 
coupling according to three (not mutually independent) categories: 
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(i) The truncation of the perturbation theory series in a s . This is sometimes considered to 
be the largest uncertainty, and it is usually estimated by applying different resummation 
techniques, such as fixed-order perturbation theory (FOPT), contour-improved fixed- 
order perturbation theory (CIPT) or renormalon chain resummation (RCPT). This has 
been studied, e.g., in Refs. El El US QUI 03] . 

(ii) The truncation of the Operator Product Expansion (OPE). The OPE is an expansion 
in inverse powers of the momentum, each multiplied by a vacuum condensate of in- 
creasingly higher dimension, which in any practical application has to be truncated at 
a certain order [12\ I13j . Furthermore, the divergence of the perturbative series in a s 
may be linked to the OPE through renormalons [Hj. This suggests a strong correlation 
between perturbation theory and the OPE. In particular, the gluon condensate depends 
on the renormalization scheme and may also depend on the order at which the series in 
a s is truncated |15| . However, if the same prescription for the perturbative expansion 
is used in both vector and axial-vector channels, the gluon condensate should also be 
the same in these two channels. The analysis of Refs. [21 [3] currently finds two different 
values for the gluon condensate from the vector and the axial-vector channels, given, 
e.g. in CIPT, by 



which are incompatible within about 3 standard deviations. This suggests that the 
systematic errors may not be fully understood. The situation concerning condensates 
of higher order is also unclear. See for example the recent analyses in Refs. [T6 | [T71 [T8l 
H51I2D1IZII 1521125], and references therein. 

(iii) The OPE is valid in the euclidean region, where there are no physical states. In contrast, 
on the Minkowski axis the expansion is not guaranteed to work. This lack of convergence 
of the OPE on the Minkwoski axis is usually referred to as due to "Duality Violations" 
(DVs). To our knowledge, the best way to understand the problem is at N c = oo, 
where the spectrum becomes an infinite set of poles. Even in this limit the OPE has 
a cut due to the logarithms from the anomalous dimensions that survive in the large- 
N c limit, and does not reproduce the set of infinite poles. Nevertheless, the standard 
analysis of tau decays involves an integration over a circle in the complex momentum 
plane which includes the Minkowski axis as well (see Eq. (jl.8p below). Of course, it is 
not known how much things might improve in the real-world case, in which N c = 3, 
relative to the limit N c —* oo, because at N c = 3 the infinite set of poles becomes a 
cut. However, it is unlikely that the OPE will reproduce this cut systematically. Of 
course, this problem has been widely recognized, and strategies have been introduced 
to ameliorate the problem. One of those strategies employed in both the ALEPH and 
OPAL analyses is the use of so-called "pinched weights," in which zeroes in the weight 
function suppress the region near the Minkowski axis [23] (see below). While this is 
expected to help, it is unknown to what degree, because there is no systematic theory 
of DVs [25]. This makes it very difficult to estimate the associated systematic error. 

In the method of analysis that has become standard for tau decays DVs are set to zero by 
fiat pQ. We consider this situation to be unsatisfactory; if eventually DVs are to be dismissed, 
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Figure 1: Analytical structure of n(q 2 ) on the complex q 2 plane. The solid curve with an 
arrow shows the contour to which Cauchy's theorem is applied. 

it should be as a result of a data-driven analysis rather than as a result of not knowing how 
to take them into account. Our aim, then, is to get a better understanding of the errors 
associated with DVs, as well as their interplay with the use of the OPE. Lacking a theory of 
DVs that one could use to perform a systematic study in tau decays, the best one can do is 
to resort to a model, and that is what we will do in this paper. 

The advantage of a model is that one has full analytical control. Therefore, it can be 
used to test the degree of accuracy obtained with the standard method of analysis for tau 
decays. The disadvantage is, of course, that the model may not represent the situation in 
QCD. Although the usefulness of a model is thus limited, we believe that a physically well 
motivated model can still give a fair idea of what to expect in the real world. Compared to 
the standard analysis, in which DVs are not taken into account at all, our model is in fact 
rather conservative, with DVs which are exponentially suppressed at higher energies. Other 
reasonable models exist in which DVs are only power suppressed and which therefore may 
well produce larger effects [25] . 

Let us start setting up the framework for our discussion. The main objects of study in 
hadronic tau decays are the reconstructed spectral functions for the vector and axial hadronic 
channels. The vector and axial current two-point correlators Hv,a(q 2 ) are analytic everywhere 
in the complex q 2 plane except on the physical axis. Use of Cauchy's theorem yields (see 
Fig. [I]) [26] 

/ ° ds P(s) -Imn(s) = — / dq 2 P(q 2 )n(q 2 ) , (1.7) 

Jo 7T 2m J\ q 2\ =S0 

where P(q 2 ) is an arbitrary analytic function, and II is either Ily or II^. As it stands, 
Eq. (jl.7p is, of course, exact. However, if the radius so is large enough, one may consider 
using the approximation II ~ H OPE on the circle part of the contour. Since, as we argued 
in point (hi) above, the OPE does not systematically reproduce the exact function over the 
full circle, a correction term appears. We will denote this correction term by T>^ p \ resulting 
in the equation 

/ ds P(s) -ImU(s) = -— & dq 2 P(q 2 )U OPE (q 2 ) + V^(s ) . (1.8) 

JO K 2m J\q 2 \=s 
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The extra term T>^ p \sq) represents the presence of duality violations. In all analyses of r 
decay data to date, "D^i (so) has been assumed to vanish. In particular, the systematic error 
reported in Eq. (jl.2p only includes the convergence of perturbation theory (point (i) above). 

One way to try to reduce the impact of non-perturbative effects is by going to the highest 
possible energy available, i.e., so = m^. Following this strategy, first put forward in Ref. [24"] . 
Refs. [2j [6] consider a family of spectral moments of the vector and axial-vector correlators 
of the form 

Rki oc / ds w k i{s) -ImII(s) , (1.9) 
Jo * 

where the "pinched weights" Wki(s) are polynomials of the form 




(1.10) 



The first factor vanishes for s = and therefore suppresses the contribution of duality 
violations near the physical axis, where their impact is expected to be the largest. The second 
factor places more weight on regions which are closer to the tau mass with increasing I, where 
perturbation theory should work best. The third factor, as well as two powers of (1 — s/m*) 
come from kinematics. A convenient set of weights consists of the five combinations wqq, W\q, 
it'll) if 12 and u>i3, which are used to perform a combined fit to extract the four unknowns 
given by a s and the condensates of dimension D = 4, 6 and 8. Although these weights are 
sensitive to condensates up to dimension 16, both the DV term T>^ kl \so) in Eq. (|1.8p as well 
as the condensates of dimension higher than 8, have to be neglected in the fit in order to 
have less parameters than data points. 

As we will see, in spite of questions about higher-order terms in the OPE and the possible 
presence of DVs, the model we will present in Sec. [3] shows that the results one obtains 
with the method of pinched weights do fairly well in determining the central values for the 
coefficients of the OPE, including the coefficient of the identity operator, which is the one 
used in the QCD analysis for extracting the value of a s . However, this does not come about 
without an associated systematic error. In the model, this systematic error turns out to 
be larger than the final error quoted in Eq. (jl.2|) . suggesting the possibility that the latter 
may have been underestimated. Moreover, if the model does not grossly misrepresent the 
real case, it also shows that, in order to make a clear improvement in the reduction of this 
systematic error, at least up to dimension-10 condensates as well as duality violations have 
to be taken into account. 

This paper is organized as follows. In Sec.[2jwe summarize briefly the theoretical analysis 
of hadronic tau decays. Section[3]is devoted to introducing the model for the vector correlator. 
We work out its high-energy expansion, and determine the duality violating terms analyt- 
ically. In Sec. H] we investigate the commonly-used "moments method" within our model. 
In Sec. [5] we propose an alternative method for the determination of the coefficients of the 
OPE, which allows us to take higher-dimension operators and DVs into account, and we test 
this on our model, comparing with the moments method. We focus on the determination of 
the coefficient of the unit operator because this is the one used for the extraction of a s in 
the real-world case. Our conclusions and an outlook for future directions are summarized in 
Seel 
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2 Summary of theoretical framework for hadronic tau decays 

Extraction of experimental information on the inclusive hadronic tau decays to vector and 
axial-vector channels proceeds through the analysis of the weak non-strange hadronic decays 
r~ — > {V~ / A~)v T . It is convenient to work with the quantity 

n^c^w (21) 

in which the hadronic decay widths are normalized to electron decay. This quantity can be 
expressed as [1] 

= 12,S EW \V ud \ 2 j™ % (l - JL) 2 [(l + 2-^) ImnW (s) + Imnj?' (.) 

(2.2) 

where H-y\ are the longitudinal (J = 0) and transverse (J = 1) components of the vector 
and axial-vector vacuum polarization functions, defined by 

njj^fo) = i /dW^(0|T{J M (*)jt(0)}|0) 



= (q»q» - q 2 9 n U^ A (q 2 ) + q»q» U { ^ A (q 2 ) , (2.3) 

in which JY(x) = u(x)^ ^(x) and J^(x) = u{x)^ ^^{x) . 

For simplicity, we will consider only the vector correlator, in the isospin limit. This means 
that we set Imlly^(s) = 0, rendering the vector correlator purely transverse. Furthermore, 
we will also take the factor S^wl K«i| 2 = 1 since this factor can always be reinstated if needed. 
With these changes, we define the moments Rki(m 2 ) as 

Rki(m 2 T ) = Utt 2 / w kl {s) -1mIL v (s) , (2.4) 

Jo m r 7T 

which can also be re-expressed, with the help of Eq. (|1.8p . as 

R kl (m 2 T ) = 6ttz / %w kl {s) Il° PE (s) + \2-K 2 V^\m 2 T ) , (2.5) 

J\s\=m? T mi 

with Wki(s) given in Eq. (jl.lOj) . In particular, R Q0 (m 2 ) = RX in Eq. (12.21) . 



3 The model 



The model we will use to illustrate the possible role of DVs follows a suggestion of Refs. 
Wt\ 128] . It is a model for the vector two-point correlator which contains an infinite set of 
resonances with a finite width. Since DVs are a consequence of the analyticity properties of 
the function Hy(g 2 ) in the complex q 2 plane, due care must be exercised to introduce this 
width in a manner that will not spoil these properties^ Specifically, the model is given by 
the expression 



2F 2 

z + m 2 p 



2F 2 



z + my(n) 



+ constant 



(3.1) 



In this regard, we note that a naive Breit-Wigner expression is not good enough because it contains an 
imaginary part even for q 2 < [25] . 
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where 

* = A 2 (^_Zi£) C , C=l-^- , m 2 v (n)=m 2 + nA 2 . (3.2) 

This model describes a separate rho meson together with an infinite tower of vector resonances 
with masses on a Regge trajectory (a "daughter trajectory"). Decay constants are given by 

{0\V lx (0)\ Pn (p,X)) = V2F mv(n)e^(p,X) , (3.3) 

and widths by 

T v (n) = ^m v {n) + o(j^j . (3.4) 

Similar expressions hold for the rho meson. For illustration, factors of N c are inserted into 
these expressions such as to allow for a comparison with the large N c limit of QCD. However, 
in our analysis, we will take iV c = 3. For a detailed discussion of properties of this type of 
model, we refer to Refs. [29], [30] . 

The function (13. lj) has a cut along the positive real axis in the q 2 plane and satisfies the 
once-subtracted dispersion relation 

r°° 1 

IV(g 2 ) = q 2 / ^—-lmU v (t) + H y (0) . (3.5) 

The (infinite) subtraction constant Ily(O) has no imaginary part and will play no role in our 
discussion since it does not contribute to equations such as fjl . 71) . We will neglect it in the 
following. 

Using the relation 

N ( 1 1\ 1 

lim V — = -Hz) 1E (3.6) 

TV^oo \ n + z n I z 

n=l x y 

where ip(z) = T'(z)/T(z) is the Euler digamma function and 'Je is the Euler-Mascheroni 
constant, it then follows that the full correlator can be expressed in closed form: 



ny(g 2 ) = i 



ip I — t-= — + constant 



z + m 2 A 2 V A 2 



(3.7) 



Figure [2] shows a comparison of the data from the ALEPH collaboration [31J and the 
imaginary part of Hy^ 2 ) in Eq. (|3.7|) with the following choice of values for the parameters 
of the model: 

F p = 138.9 MeV , m p = 775.1 MeV , F = 133.8 MeV , = 0.158 . (3.8) 

With these values, Eqs. (|3.14|3.15p below determine A and mo to be 

A = 1.189 GeV , m = 1.493 GeV . (3.9) 

Although the model is clearly far from perfect, Fig. [5] shows that it also does not grossly 
misrepresent the main features of the spectrum, even though the value for the parameter £ is 
very close to unity (£ ~ 0.95). Therefore, the model illustrates how the large N c expansion 
can be quite consistent with a realistic description of the real-world spectrum, for which 
A^ c = 3. 
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Figure 2: Experimental data for the spectral function for the vector channel from the ALEPH 
collaboration j2[ [31]. In blue, we show the curve given by the model defined in Eq. (|3.ip with 
values for the parameters as given by Eq. (|3.8j) . 



3.1 Operator product expansion 

Using the asymptotic expansion of the digamma function 



1 °° B 

i>(z) =logz- — , \ z \ » 1 , -7T < arg(z) < vr , 

71=1 



(3.10) 



where B2 n stand for the Bernoulli numbers Bq = 1,B\ = — \,Bi = ^, etc., it is straightfor- 
ward to obtain the OPE for the correlator (13.11) as 



^p- Co log(-<? 2 ) + ^2 ~T~ + constant 

k=l Z 



(3.11) 



where [29 



C n = l 



C 2fe = - ( (-l) k 



/ ;; rn^ + ^A^ F z B k l 



(3.12) 



The latter play the role of the vacuum condensates. In the previous expression, B^ (x) stand 
for the Bernoulli polynomials 



k=Q ^ ' 



(3.13) 



Comparing with the case of QCD, we see that if we make the identification 

El 

A 2 24^ 2 ' 



(3.14) 
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the coefficient of the parton model is reproduced. Furthermore, imposing the constraint 



F: 



A 2 



(3.15) 



we require that the condensate of dimension two vanishes, i.e., C2 = oH 



Let us now look at the contribution from the OPE to the moments. Defining Rffl as the 



contribution of the dimension- D operator to the moment R^i, 
a straightforward calculation yields (k>l,j>0): 



R 



(0) 

00 

(0) 

1? 



R 



(2k) 
00 

(2k) 



12tt 2 F 2 

A 2 
144vr 2 F 2 

A 2 
72tt 



Cn 



Cn 



i.e., Eq. ([23]) with T>W = 0, 



(3.16) 



7 + 3j 



(m 2 ) k< > 
72vr 



a 



(5 + i)(4 + j)(3 + i)(2 + j)(l+j) 
A 2fc(C-l) sin^vr^ 



2/v 



1/ 



(fcC-4)(fcC-3)(fcC-l) ' 

A 2fc(C-l) ( 3/tC _ 



(m 2 ) fc C ~ 2fc (k( - 5 - - 4 - - 3 - - 2 - j)(k( - 1 - j) 



7 — 3j) smkir( 



In the limit N c — * 00 (i.e., C — > 1), the contribution of the condensates to a given moment 
is only non-vanishing whenever the dimension 2k is such that there is a zero in k in the 
denominator of the corresponding expression for the moment. The moments, however, are 
well defined for all values of k. For instance, for Rqq'^, the only non- vanishing values happen 
for k = 1, 3 and 4, which correspond to dimensions D = 2, 6 and 8, in the limit £ — > 1. Since 
in the model £ is very close to one, these terms also dominate in the case of finite A^ c . 



3.2 Duality violations 

Up to now duality violations have not made an explicit appearance. However, the restriction 
on the validity of the asymptotic expansion (|3.10p to the region | arg(z)| < tt implies that the 
OPE (13. lip cannot be valid on the full contour \q 2 \ = sq, causing a breakdown of the naive 
application of Cauchy's theorem in which one would freely replace H(q 2 ) by IL OPE (q 2 ) in 
Eq. (11. 7p . The contribution from the Minkowski axis (arg(z) = n <-> Reg 2 > 0) is responsible 
for this breakdown. This is how duality violations make their appearance and why T>^ kl \m 2 ) 
in Eq. (12. 5ft does not vanish in our model. 
Use of the reflection property 



lp(z) = 1p( — z) — 7TCOt (m) 



(3.17) 



allows us to write an expansion for large values of \q 2 \ which is valid for Reg 2 > (in 
particular on the Minkowski axis): 



with the correction A given by 

2ttF 2 1 



U OPE^2 



A(q 2 ) 



A 2 C 



-i + COt 7T 



A 2 



+ 7T 



A 2 



(3.18) 



(3.19) 



2 In QCD this would mean that we neglect terms of order m 2 q /z relative to those of order m q (tptp)/ 



s 



This function A(g 2 ) behaves for large values of complex momentum q 2 as |29| 



A(, 2 



~ e 



27r I "ia{(T-*)C}l 



€ [0,7r/2] U [37r/2,27r) . 



(3.20) 



This behavior shows that the limit N c — > oo (i.e., £ — > 1) and the limit (ft —* do not 
commute. If the limit — ► oo is taken first, the expression (|3.19p reveals the existence of 
poles at <j) = every time that q 2 = m 2 , + nA 2 . This is as it should be, since in this case the 



spectrum of the model (13. ID is composed of an infinite set of Dirac deltas located precisely 
at these points. For finite N c , on the Minkowski axis, one obtains instead 



cos 



- ImA(t) 

7T 



27T 



t \< 



COs(7rC) + 



-27rsm[7rC](^-J 



C A2 cosh 2tt (^) f sin(vrC) -cos 2vr ((^) C cos(vrC) + $ 



A F Mil 

CA 2 



sin(7r£) 



cos 



2^1 (^) C cos«) + ^ 



(3.21) 



(3.22) 



where the second expression is valid for t large. 

The DV function T>^ kl \ which corrects for the naive use of the OPE in Cauchy's theorem 
as in (|2.5p . is given by (see Appendix and Ref. [29] ) 



2>M, 



OO j 

ds w k i(s) — ImA(s) + O | r 

,2 7T 



2tt- 



(3.23) 



It follows from Eq. (pT22j) that the DV function V^ kl \m 2 T ) behaves like ~ e 2?r ^ ^ . There- 
fore, duality violations are exponentially suppressed in our model, and thus parametrically 
smaller than any term in the OPE. However, we also note that the exponential suppression is 
not that fast, because of the presence of the factor a/N c « 0.158 for the choice (|3.8p . Thus, 
DVs may be numerically important at or below the tau mass, and the question is, by how 
much? 

In order to answer this question, we would like to study to what extent one can determine 
the OPE coefficients, given a certain analysis method. To be more specific, consider our 
model, cf. Eq. (13. lip . If one starts with a general parametrization of the OPE of the form 



n 



OPE-m 

V 



(A2fc;<? 2 



2F 2 
A 2 



A log(-<? 2 



C2k 

+ Y1 X2k ~k ' 

k=i z 



(3.24) 



of course, the exact answer is given by Tly PE ~ &t (\2 k ,q 2 ) with Ao = A2 = A4 = . . . = 1. How 
many of the parameters \2k,k = 0,1,2,... can one determine, and with what precision^ 
What is a good method to minimize the error | A2A; — 1| in this determination? 



4 The moments method applied to the model 

We are now in the position to test the standard method for determining the OPE coefficients, 
i.e., the X's in Eq. (pOIjh 

3 Sincc in our model we set the dimension-2 condensate C2 = 0, this means that the parameter A2 is 
irrelevant. We will disregard A2 in what follows. 
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OPE,ji m _ 


_ 8 OPE dim _ 8 (DV) 


IFdim-8 


IF(jim-10 


IF di m-io(no DV) 


|A -1| 


0.016 


0.017 


0.025 


0.00077 


0.0095 


A 4 


0.85 


0.85 


0.56 


0.96 


1.08 


A 6 


0.65 


0.65 


0.56 


0.94 


1.05 


A 8 


0.53 


0.54 


0.31 


0.84 


0.97 


Aio 








0.59 


0.72 




Table 1: 


Results of various 


fits discussed in Sees. 


[1 and El 



As mentioned in the introduction, this method consists of a global fit to the spectral 
moments R T = i?oo> -RiOi -Riii R12 and R13, defined in Eq. (|2.4|) . evaluated at the tau mass, 
and taking into account only the contributions from the OPE up to D = 8. Of course, in the 
real world the parameter Ao is actually an expansion in powers of a s of the form 

A = l + — + .... (4.1) 

7T 

Since our model is too simple to have any perturbative corrections, we will take the value 
obtained for |Ao — 1| as a rough estimate of the accuracy reached in determining a s /7r. 

Following Refs. [21 EH IE]; we compute the five spectral moments at the tau mass as a func- 
tion of the parameters A 2 fc with k = 0,2,3,4, using Eq. (|3.24|) . and making the replacement 
C 2 k -> A 2 fcC 2 fc in Eqs. (13T61) . The DVs, i.e., the functions V^- kl \m 2 T ) in Eq. ([23]), as well 
as A 2 fc for k > 4, are set equal to zero. We then extract the values of A 2 fc,/c = 0,2,3,4 by 
doing a combined least-squares fit to the five moments obtained from the "data" in Eq. (|2.4p , 
in which the true spectral function ^Imlly(q 2 ) of Eq. (|3.7p . with parameters chosen as in 
Eqs. (|3,8|3.9p . is used. Although within our model we have access to an infinite number of 
points, we have discretized the spectral function for the model in the same way as the ALEPH 
data have been binned, i.e., into 125 bins with a binsize of 0.025 GeV 2 each, with starting 
and ending points at sq = 0.0125 GeV 2 and so = 3.1375 GeV 2 , respectively. 

The results for the fit are shown in the second column of Table[TJ As one can see, the values 
obtained for the A's are fairly close to unity (in particular for the "perturbative expansion" 
Ao), showing that the method works reasonably well. However, if the analogy (|4.ip is used, 
this translates into a systematic error in a s of 15%, i.e., larger than the estimated theoretical 
error shown in Eq. (jl.2p . Of course, this does not mean that the analysis of Refs. (2j 3, 6] must 
have a systematic error of this size. However, it does suggest that this method of analysis is 
not safe against non-perturbative effects beyond those included in the ALEPH analysis. 

Of course, in the context of our model, part of the reason may be that DVs have not 
been included, since the functions p[ fc ^(m 2 ) are actually not zero. This is why we have 
repeated the analysis with the same five moments and condensates, but now including the 
exact duality violations of the model, i.e., Eq. (|3.23[) calculated with the exact expression 
(|3.2ip . In other words, we first "subtracted" the exact DVs from the "data" generated from 
the exact spectral function before extracting the parameters A 2 fc. The results are listed in 
the third column of Table [1] As one can see, inclusion of duality violations in the analysis of 
moments at m T does not by itself lead to an appreciable change in the results. Consequently, 
we conclude that also the effect from higher dimension condensates may have to be taken 
into account. 
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In fact, as already mentioned in the Introduction, the moment -R13 is sensitive to con- 
densates with dimension up to 16. The structure of the pinched weights (ll.lOp . while being 
good at canceling the integration region near the Minkowski axis, also makes the moments 
sensitive to condensates of higher order, especially for larger I. However, if we only use the 
moments at the fixed scale mif, there is simply no way to include more condensates while 
using the same moments and still have a "fit" in which there are more data points than 
degrees of freedom. 

5 An iterative fitting method 

The reason for using the values of the moments at the scale of the tau mass is that non- 
perturbative effects are smaller at a higher scale, and that includes DVs as well. However, 
in order to include the DV effects, we will have to go to lower scales to be able to have more 
data points, as we will explain below. To this end, let us imagine that the tau mass is an 
arbitrary scale we could change at our convenience, and let us refer to this scale as sq. The 
moments are now given by the same expressions (|2.4|) and (|2.5p . but with the replacement 
m\ — > sq everywhereEl This trick has also been used by ALEPH and OPAL when studying 
the validity of perturbation theory as a function of the scale so [2J [6] . 

Based on experience with applying our model to Hy — n^, we propose a simple parametriza- 
tion of the DVs through the ansatz |29j 

-Im A Dy - fit (s) = k e~ 7S sin (a + /3s) , (5.1) 

in which a, (3, k and 7 are parameters to be fitted. This ansatz is of course inspired by the 
asymptotic behavior of ImA given in Eq. (I3.22p pl We should thus be able to fit it to our 
model to see if the results for the coefficients of the OPE improve relative to the standard 
analysis described in the previous section. 

Naively, one would think of fitting the moments to extract the values for all the OPE 
coefficients X 2 k together with the DV parameters a, /3, 7 and k. However, this leads to a 
complicated, multi-dimensional fit which yields poor and unstable results for the parameters. 
Furthermore, the pinched weight moments have been designed not to be very sensitive to 
DVs, and thus the determination of the DV parameters from the moments is difficult. On the 
other hand, the spectral function is sensitive to DVs but less sensitive to the OPE coefficients 
because they only enter through logarithms. Except for the parton model logarithm, these 
logarithms are always screened: in the real world by an extra power of a s and in the model 
by a factor 1 — C- Consequently, it appears to be best to combine using the spectral function 
for fitting the DVs with using the moments for fitting the OPE coefficients. 

Therefore, the strategy we propose consists of an iterative two-step procedure. Our 
starting ansatz for the spectral function is 

-lmU^(X 2k ,s) = ±lmn° PE - Rt (\ 2k ;s) + - Im A Dy ~ fit ( S ) , (5.2) 

7T 71" 7T 

with the starting values X 2k = 0, k > 1, and the functions Hy PE ~ {it (X 2 k; s) and ImA Dl/_fit (s) 
given by (|3.24p and (|5.1|) . respectively. Note that Ao is not set to zero in Eq. (j 5 . 2 j) to account 
for the parton model contribution, which is the dominant part of the OPE and not suppressed 

4 This includes replacing — > so in Eq. (|1.10|l . 

5 In the real QCD case, the function (|5.ip will have to be compared to the experimental data. 
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•Smin (GeV 2 ) 


least-squares 


|A Q - 1| 


A 4 


A 6 


As 


Aio 


1.0125 


7.4E-5 


0.01181 


0.828 


0.808 


0.701 


0.472 


1.2125 


1.8E-5 


0.00589 


0.901 


0.876 


0.778 


0.531 


1.4125 


8.2E-7 


0.00077 


0.963 


0.935 


0.843 


0.590 


1.6125 


3.3E-6 


0.00026 


0.975 


0.947 


0.857 


0.602 


1.8125 


2.5E-6 


0.00223 


0.934 


0.903 


0.794 


0.529 


2.0125 


2.6E-5 


0.00246 


0.920 


0.886 


0.762 


0.485 



Table 2: Results of the iterative fit as a function of s m i n ; see text. 



by a s in the real world (or by 1 — £ in our model). As before, we discretized the spectral 
function for the model into 125 bins with a binsize of 0.025 GeV 2 each, with starting and 
ending points at so = 0.0125 GeV 2 and so = 3.1375 GeV 2 , respectively. 

In a first step, the ansatz (|5.2j) is fitted to the "data" given by the exact spectral function 
ilmlly(s) over the interval 1.4125 GeV 2 < s < 3.1375 GeV 2 . Values for the parameters 
Ao, a, (3, k and 7 are obtained from this least-squares fit. Once these values are known, the DV 
function T>^ kl \so) can be calculated from Eq. (I3.23h with Eq. (15. lj) . and then be incorporated 
into the moments Rki(so) of Eq. (|2.5|) . 

In a second step, the values Ao,4,6,8,io & r e then determined from a second least-squares fit 
to the Rki(so) moments, which include the duality violations determined in the first step, as 
a function of the scale so in the window 2.0125 GeV 2 < sq < 3.1375 GeV 2 . Note that a, f3, 
7 and k are held fixed in this step, while Ao is of course allowed to vary. The values obtained 
for the OPE parameters A2/C, for k = 0, 2, 3, 4, 5, are then held fixed while the function (|5.2p is 
re-fit again in order to find a new set of values for a, (5, 7 and k. This second step is iterated 
until stability is found. The choice of the lower end in the above fitting windows is dictated 
by the stability of the result as well as by the minimum value of the least-squares function. 
The upper end is fixed by the tau mass, as it would be in the real QCD case. 

Table [2] shows the result of the fit to the least-squares function constructed with the 
three moments Rqq, -Rio an d R\\ as a function of the lower end in the fit of the function 
(|5.2p to the spectrum, together with the values of Ao,4,6,8,io- As one can see, the smallest 
value of the least-squares function is obtained for s m i n = 1.4125 GeV 2 , which we therefore 
take to be the "best" determination. The values of Ao,4,6,8,io f° r this fit are also shown in 
the (bold-faced) fifth column of Table [TJ Note that the values at s m i n = 1.6125 GeV 2 are 
actually closer to the exact values in the model, but the least-squares function is larger than 
at s m i n = 1.4125 GeV 2 . The spread of the values in the parameters Ao,4,6,8,io between these 
two cases can be taken as a rough indication of the systematic error of our iterative method. 

For the fit in the second step, we have tested different combinations of spectral moments. 
The best determination in terms of stability and the smallest value of the least-squares 
function occur for the fit with the three moments Roo, Rio and R\\. The fits gradually 
worsen if one uses more or less moments. The lack of improvement when adding extra 
moments agrees with the fact that higher-order moments are sensitive to higher orders in the 
OPE. Therefore, the inclusion of higher moments is disfavored because of the OPE truncation, 
as we already anticipated. The lack of improvement with too few moments is related to the 
fact that in that case the OPE parameters are too correlated and not sufficiently constrained 
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Figure 3: Comparison between the exact spectral function (black solid curve) and the iterative 
fit described in Sec. [5] (blue dashed curve). 



to be determined by the fit. Therefore, the three-moment fit happens to be the optimal choice 
in the model. 

Figure [3] shows the final estimate for the spectral function Eq. (|5,2p obtained with the 
parameter values from our iterative fit (blue dashed curve), in comparison with the exact 
spectral function (black solid curve). We note that the rho is far away from the iterative 
fit. This happens because the rho is not part of the tower of resonances described by the 
digamma function in Eq. (|3.7p ; in other words, the rho is not part of the Regge trajectory. 
This does not cause a problem, since the approximation (15. ip is only needed for s > sq, 
cf. Eq. (|3.23p . where we always take so > 2 GeV 2 . The final values of the DV parameters 
extracted from the fit are 

k = 0.035, 7 = 0.51 GeV -2 , a = 12, (3= -4.3 GeV" 2 . (5.3) 

These numbers may be compared to the results obtained in the model, if we make the 
approximation « t in Eq. (|3.22p . Since £ ~ 0.95, this approximation is justified. One finds 

K = 0.051, 7 = 0.71 GeV" 2 , a = 10, = -4.5 GeV" 2 , (5.4) 

showing that the iterative fitting procedure does quite well for our model. 

As already mentioned above, for the determination of the OPE coefficients, i.e., the 
parameters A2fc, the result of the fit is summarized in the fifth column of Table [TJ For 
comparison, we also give in the fourth column the result we obtained if the condensates 
of dimension 10 and higher are not included. Finally, the last column shows the results 
when dimension-10 condensates are included, while at the same time duality violations are 
neglected (i.e., with k set equal to zero). In this case, we note that the higher condensates 
actually move somewhat closer to their exact values. However, the value obtained for Ao, 
which is the parameter most directly related to a s , cf. Eq. (|4.ip . is substantially worse. 

Based on this table we draw the following conclusions, valid for our model. The inclusion 
of duality violations is necessary in order to make a clear improvement in the determination of 
the different terms in the OPE. However, this is only true provided enough of these OPE terms 
are included in the analysis as well. In the case of the model, this requires the inclusion of up 
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to dimension-10 condensates. We have checked that including dimension-12 condensates in 
the analysis does not improve the determination of the parameters \2k- However, the optimal 
number of terms in the OPE may depend on the case at hand, and does not have to be the 
same in the case of QCD. 

6 Conclusions and outlook 

The present analysis of data for hadronic tau decays is based on the use of finite-energy sum 
rules, in order to deal with the non-perturbative effects that cannot be neglected at a scale 
as low as the tau mass. This involves an extrapolation of the OPE for the hadronic spectral 
functions over a full circle around the origin in the complex momentum plane, including a 
region near the Minkowski axis where the physical spectrum is located and the OPE is not 
expected to work. Although the contribution from this region is minimized by including in the 
integrand convenient weight functions which vanish where the circle intersects the Minkowski 
axis, contributions from duality violations do not exactly vanish, and it is reasonable to expect 
that they have a residual effect. This effect, if sizeable, contaminates the determination of 
a s and the OPE condensates from hadronic tau decays. 

In fact, there already seems to be some (perhaps circumstantial) evidence that there is 
something beyond the OPE lurking in the tau data. As emphasized in the Introduction, 
the value of a s extracted from tau decay is barely consistent with the other most precise 
determination, which is claimed by the lattice. Moreover, the gluon condensate, which is the 
leading non-perturbative effect in the OPE, comes out clearly different depending on whether 
the vector or the axial- vector channels are used in the analysis. 

In spite of this, to the best of our knowledge, no analysis of duality violations has been 
done on the tau data, and this gap has been the main motivation for the present work. With 
the help of a physically motivated model, we have studied a method for extracting both the 
coefficient of the parton model logarithm (from which a s is extracted in the real case as a 
radiative correction), as well as the OPE condensates. For simplicity, we have taken the 
vector channel as the case study, but the axial-vector channel can in principle be treated 
similarly. The obvious advantage of having a well-defined model is that the exact answer is 
known. Nevertheless, we have tried to simulate a semi-realistic situation by imagining fitting 
a general parametrization of duality violations to the model, as if it provided the real "data." 
This allowed us to first test the standard method that has been used to date for analyzing 
tau decay data. We have found that, taken at face value, this method does a fair job, in spite 
of the presence of duality violations. The associated systematic error, however, would be too 
large to be acceptable in the case of the real world, if the results obtained with our model 
are any indication. Needless to say, we do not yet know if that is the case. 

It is interesting to note that our model provides a possible explanation for the fact that 
even if duality violations are exponentially suppressed, they might still be numerically im- 
portant at energy scales of order a few GeV. As discussed in Sec. 13.21 the scale for this 
exponential decay involves a factor N c , and may thus differ substantially from the typical 
hadronic scale of 1 GeV. 

In order to improve on the systematic error, we have proposed a new method of analysis 
based on an iterative fit, described in Sec.EJ Our results show that a substantial improvement 
may be possible if one is willing to do some modeling of the duality violations in the form of 
Eq. (|5.ip . At the same time, our results indicate that also enough terms of the OPE should 
be taken into account. In the model, we needed to include condensates up to and including 
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dimension ten. In the real world, this number of course could be different. 

Since there exists at present no systematic theory of duality violations, we had to resort to 
a model and one might thus dismiss the results on the basis that they are model dependent. 
However, we think that it is very important to realize that, by the same argument, all the 
analyses of tau data performed so far have also been done with a model, namely one in which 
duality violations are dismissed by fiat, and not by a fit to the data. In terms of our ansatz, 
Eq. (|5.ip . current analyses are equivalent to setting the value of k equal to zero. 

Of course, it may turn out to be difficult to apply our new iterative fitting procedure 
to real-world data, because these are afflicted with sizable and correlated errors which we 
completely ignored in the present paper. In addition, in a fit over an interval such as proposed 
in Sec.[5]applied to real data, the running of a s would have to be taken into account. However, 
again, that does not justify ignoring duality violations altogether. It would be interesting to 
see how our method fares in the case of real data, even if it can only be used to establish 
bounds on the parameters k and 7. 

Note added: While this work was being finalized, two new determinations of a s at the tau 
mass appeared |32t I33j. which, after running to the Mz scale, translate into 



respectively. These results differ from the lattice result (c/. Eq. (jl.5p ) by 1.5 and 2.6 standard 
deviations, respectively. 
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In Sec. [2] we already mentioned that the spectral moments can be related to the large-g 2 
expansion of the vector correlator as 



a s (M|) = 0.1202 ± 0.0019 , 
a s (M 2 z ) = 0.1212 ±0.0011 , 



(6.1) 
(6.2) 



A Derivation of Eq. (p3j) 



Rkl{m 2 T ) = 6ni 



ds 



1 w kl (s)U PE(s) + l27r 2 V^(m 2 T ) . 



(A.l) 



\so\=rri: 



2 m~ 
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Figure 4: The contour used to evaluate Eq. (IA.2j) . see text. 

The OPE part was discussed in Sec. (|3.ip . The second term on the right-hand side corrects 
for the fact that the OPE is not valid on the physical axis, and it can be expressed as 

rt kl \s,)=! ^w kl (q 2 )A(q 2 ), (A.2) 

J\q 2 \=s , Req 2 >0 m r 

where the contour runs counterclockwise only over the right half hemisphere, cf. the bold- 
faced semi-circle in Fig. H The function A(q 2 ) is defined in Eq. (ETTUl) . 

Unlike in QCD, the model allows us to calculate the duality-violating part as the difference 
between the full correlator and its analytically-continued OPE expansion. Using the reflection 
property of the Digamma function (|3.17p . one easily obtains Eq. (|3. 19j) . In particular, note 
that the resonance poles appear in A(q 2 ), as can be seen through the identity 

7r cot (irz) = - + 2z — (A.3) 



^-j z 1 

n=l 



Expanding the cotangent for large q 2 and small a/N c , one obtains Eq. (|3.2U|) . Note that for 
small a/N c , the damping factor along the imaginary axis ((f) = tt/2,3tt/2) is stronger than 
along the real axis ((f) = 0): 



7T 37t\ _ 2 7T 



o \q 2 \ a 

A((f> = 0) ~ eT^-^m . (A.4) 

Due to this difference, the leading contribution to the spectral moments can be easily derived 
by closing the contour as shown in Fig. 2] and using Cauchy's theorem. Because of the 
exponential fall-off, the integral vanishes at infinity, and we are left with integrals over the 
real and imaginary axes: 

T>W(a ) = - dsw kl (s)-ImA(s + ie) + \ -/ \ dq 2 u; kl (q 2 ) A(q 2 ) . (A.5) 

J so ^ \J — iso Jiso J 

Using the fact that the integral over the real axis dominates, cf. Eq. (|A.4|) . one finds the 
result given in Eq. (|3.23p . 
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